Quantum tomography for solid state qubits 
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We propose a method for the tomographic reconstruction of qubit states for a general class of solid state 
systems in which the Hamiltonians are represented by spin operators, e.g., with Heisenberg-, XX Z-, or XY- 
type exchange interactions. We analyze the implementation of the projective operator measurements, or spin 
measurements, on qubit states. All the qubit states for the spin Hamiltonians can be reconstructed by using 
experimental data. 
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Quantum information processing requires the effective 
measurement of quantum states. However, a single quan- 
tum measurement can only obtain partial information of a 
quantum state. The reconstruction of a quantum state re- 
quires measuring a complete set of observables on an en- 
semble of identically prepared copies of the system. This 
method called quantum state tomography [1], is very impor- 
tant because any unknown state can be ascertained by to- 
mographic measurements. Moreover, the full description of 
qubit states can increase the accuracy of quantum operations. 
Tomographic measurements have been experimentally imple- 
mented for, e.g., the nuclear spin state of an NMR system |2], 
the electromagnetic field and photon state |3|], the vibration 
state of molecules |4], the motional quantum state of a trapped 
atom 1 5], and atomic wave packets |6]. 

Experimental investigations on solid state qubits are very 
promising, especially in superconducting 0, [U and quantum 
dot structures (2J. These recent achievements make it nec- 
essary to experimentally determine quantum states in solid 
state systems. Although there are many theoretical studies 
on tomography (e.g., references 1 1 OH and references therein), 
to our knowledge, these are not specific to solid state sys- 
tems. Here, we focus on this question for quantum com- 
puting models using standard spin representations for solid 
state qubits. Our proposal is related to tomographic measure- 
ments using NMR. The measurements of the density matrix 
in NMR experiments are obtained from the NMR spectrum of 
the linear combination of "product operators", i.e. products 
of the angular momentum operators 111 ill . However, experi- 
ments in solid state systems usually involve the local single- 
qubit projective operator measurement (POM) or spin mea- 
surement. So, we study the method of tomographic recon- 
struction of solid state qubits by POM or spin measurement 
for a number of promising solid-state quantum computing 
models El El El El ■ We will investigate 
how the multi-qubit correlation measurements can be realized 
by virtue of an appropriate two-qubit operation combined with 
single qubit operations. 

State and measurements. — Using the density matrix form, 



an n-qubit state p can be expressed as 

p= Y l ^ ril "' ln ° h <g " T ' 2 (8> '" <glcr/ "' (1) 

ii,-- ,l n —0,x,y,z 

where r; 1 ...; n are 4™ real parameters, a-i m=x ,y.z and 
cz m= o (0 < m < n) are the Pauli spin and identity opera- 
tors of the rath qubit respectively. We adopt the convention 
|0) = | |) and |1) = | I) to denote the computational basis 
states of each qubit. The normalization condition Trp = 1 
makes ro,....o = 1 which means that p can be specified by 
(4" — 1) real parameters. These parameters correspond to 
the expectation values of the measurements given by the op- 
erators <Tj 1 g) • • • (g) o~ j n ; that is, TT{p(aj 1 ® •• • ® c^)} = 
ri 1 ...i n Sj 1 i 1 ■ ■ ■ where ii, • • • , l n are not simultaneously 
taken as zero. If there are n — m identity operators among 
o~i x g) • • • g) ai n , the measurement is really done by the m 
qubits and it can be abbreviated by the tensor product of 
only the m Pauli operators, which is denoted hereafter by 
o'ii 1 g) • • • g) o- m i m . The (4" — 1) measurements required 
to reconstruct the n-qubit state can be decomposed into a 
summation from the single-qubit to n-qubit measurements as 
THj=i ^ ( j) ' wnere 3 J (™) is the number of j-qubit measure- 
ments and ("J is the binomial coefficient. 

To reconstruct the n-qubit state p, we need to determine all 
of its expanded coefficients {r^ ... j n }. In solid state systems, 
the correlated multi-qubit measurement is not realizable now, 
and the experimental readout is often done via single-qubit 
POM (e.g., RefsJ0000) or single spin measurement 
(e.g., Refs. II6II17I0 . Without loss of generality, we assume 
that the POM is denoted by |1)(1| and the spin measurement 
is presented by a x or a y . Below, we will discuss how to deter- 
mine the (4" — 1) coefficients of an n-qubit state by using ex- 
perimental data of the POM and then generalize it to the spin 
case. Our goal is to build a correspondence between above the 
measurements and actual measurements done via |1) (1|. 

Single-qubit measurements. — The single-qubit measure- 
ments aiz (I = 1, 2, • • • , ?i ) can be done by the projectors 
(|l)(l|)i due to the equivalence ai z = am — 2(|l)(l|)z, with 
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identity operators er;o- Thus the POM experiments can di- 
rectly determine n coefficients via n outcomes of the mea- 
surements (|l)(l|)j. The measurements corresponding to the 
remaining (4™ — n — 1) coefficients cannot be directly per- 
formed because of limitations of current experiment. In order 
to obtain these coefficients, a sequence of quantum operations 
W is needed such that each coefficient can be transformed to a 
position that is measurable by a POM experiment. The prob- 
ability p of the kh single-qubit measurement (|1)(1|)/ on the 
operated state p can be expressed as 

V = Tr{WpW^(\l)(l\)i} = Tr{pWmi)(l\)iW}, (2) 

which means that the experimental POM (|1)(1|); on the state 
WpW^ can be equivalently described as Wt(|l)(l|)zW — 
[h ~ (Ji z W~\ on the original state p. Because the POM 
experiment is equivalent to W'ai z W, we can call it an equiv- 
alent rn-qubit measurement of (|1)(1|);, if the largest num- 
ber of Pauli operators in the operator product expansion of 
W^aizW is an m-qubit. Thus the (4" — n — 1) coefficients 
can be obtained by the single-qubit projection (|1) on the 
state p with appropriate operations W implemented by the dy- 
namical evolution of the system with experimentally control- 
lable parameters. 

For all universal quantum computing proposals, the most 
general Hamiltonian of the system can be described as 

n n 

H = ^ ^ ei a <Jl a + ^ X! JlL G la® G rali, 

1=1 a=x,y,z l=i<m a,f}=x,y,x 

(3) 

where {ei a } and { J"^ } are controllable and tunable system- 
specific one-qubit and exchange coupling parameters, which 
are required by the universality of quantum computing |20], 
07 a=x,y.z denote the Pauli operators of the kh qubit. Without 
loss of generality, all parameters are assumed to be positive 
real numbers. 

In order to obtain each coefficient corresponding to single- 
qubit measurements ai x or ai y , all single-qubit operations 
need to be performed separately by controlling the one-qubit 
parameters Ei a while turning off all interactions in the Hamil- 
tonian (0, that is, J?Z, = 0. For n single-qubit measurements 
{<Jiy}, each a\ y can be equivalently obtained by (|l)(l|)j, 
after the kh qubit is rotated tt/2 about the x axis, the lat- 
ter expressed as X[ = exp{— i7rai x /4}. This rotation can 
be realized within the evolution time t = fnr/4ei x , after the 
one-qubit parameters ei y and ei z are adjusted to zero. Other 
n single-qubit measurements {aix} can also be obtained by 
measuring (|1)(1|); on the state, within the evolution time 
t = fnr/Aeiy, after ei x and ei z are set to zero. This quan- 
tum operation is equivalent to a 7r/2 rotation of the kh qubit 
about the y axis, which is denoted by Yi = cxp{— iiraiy/A}. 
However, not all of the three one-qubit parameters Ei x , E\ v 
and £i z appear in the Hamiltonian of most physical sys- 
tems. For example, the parameter ei y is always zero in the 
charge qubit system |12|. For this case, to obtain the rota- 
tion angle 9 about the y axis we need to alternatively turn 



TABLE I: Quantum operation and two-qubit measurements. 
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on and off the single-qubit quantum operations: (1) Xi — 
exp{—isi x ai x ti/h}, with the operation time t\ = 3tnr / '4ei x ; 
(2) Zi{6) = exp{-i9a lz /2}, with 9 = e lz t 2 /h; and (3) 
Xi = cxp{~iei x (Ti x t 3 /h}, with_^3 = for/4£ fa . These can 
be expressed as Yi{9) = XiZ l (9)Xi = cxp{— i6ai y /2}. Es- 
pecially, we denote the rotation tt/2 about the z axis by Zi = 
exp{— iTTdiz/A}. In principle, if the kh qubit system has only 
two controllable one-qubit parameters e a and Ep, then the ro- 
tation angle e a /3 7 9 about the axis 7 can be obtained by first 
doing a rotation of 7r/4 about the axis a, then a rotation of 9 
about the axis (3, and finally a rotation by — n/4 about the axis 
a; that is, exp{-ie Q/37 (9cr 7 /2} = e -*W4 e -<0<rp/2 e i7r<r o /4 ( 
where a, f3, 7 can be x , y, or z and the Levi-Civita tensor e a ^ 7 
is equal to +1 and —1 for the even and odd permutation of its 
indices, respectively. To reconstruct a single-qubit state, three 
single-qubit measurements a a (a — x, y z) are sufficient to 
obtain r z , which determines the probabilities of finding and 
1, as well as r x and r y which determine the relative phase of 
|0> and 

Two-qubit measurements. — The above discussions show 
that all the single-qubit measurements can be experimentally 
implemented by POMs (|1)(1|)/ on the given state with or 
without single-qubit quantum operations. However, the im- 
plementation of multiple-qubit measurements needs non-local 
two-qubit operations. The basic two-qubit operation can be 
derived from the time-evolution operator Ui2(t) of a pair of 
coupled qubits, labelled by 1 and 2, whose Hamiltonian H12 
can be obtained from Eq. (0 with n — 2. Without loss 
of generality, we assume: (i) e\ a = £20 = £a\ and (ii) 
^im = Jhn^aP m me Hamiltonian 0, because by applying 
local unitary operations, e.g., 12111 . the Hamiltonian can al- 
ways be transformed into a diagonal form — which is actually 
used for a number of promisin g so lid-state quantum comput- 
ing models G2 H H E Hi H d . Then in the basis 
{|0i, 2 ), |0i, 1 2 ), 1 1 1 , O2), 1 1 1 , la)}, the time-evolution oper- 
ator U\2{t) is 

4 

U l2 (t) = exp{-iH 12 t/h} = Y, e ~ iEet,h \^)^9\ ( 4 ) 

9=1 

where \ip g ) (g = 1,2,3,4) are four eigenvectors of the 
Hamiltonian Hi 2 . The corresponding eigenvalues E\ = 
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— J 12 ~ J 12 ~ J 12' Ez> ^3 an d ^4 are gi ven [22] by the solu- 
tions of the cubic equation of the parameter E. Here, we only 
focus on two typical Hamiltonians which play an important 
role in the process of two-qubit operation for the most repre- 
sentative solid state quantum computing models. One is that 
all of the one-qubit parameters are switchable, for example, 
quantum dots in cavities 112 311 . However, due to technical con- 
straints and difficulties, it was found 1 24] that not all the one- 
qubit parameters are switchable in the two-qubit operation, for 
instance, for spin-coupled quantum dots 1 14], donor-atom nu- 
clear or electron spins 11811 . and quantum Hall systems 11911 . 
two one-qubit parameters such as e x and e v are switchable, 
but e z is fixed. The basic two-qubit operation with fixed e z is 



1 



U 12 (t) = -(e l0 cos 7 + e-^cos/3)/ + i 



.(l-a 2 )c _ 



i<j> 



x sin/3((7iz + a 2z ) + ~( e 10 cos /3 - e^cosj) a lz ® <j 2z 



-i-(e^ sin7 + 2ace ^ sin /3) a\ x <g> <7 2a 



*o( e ^ sin 7 — 2ace 1<p sin/3) ai y <E> &2y 



i<j> 



(5) 
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i jz 
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2b 
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1, with b 



" (J X 12 
z/ \Jl2 



Jf 2 ), and 



1/(1 + a 2 ). We also assume in Eq. that the parameters 
satisfy conditions [2 Jf 2 ± ( Jf 2 + Jf 2 )] 2 ^ 4e 2 + ( Jf 2 - Jf 2 ) 2 
and 2 Jf 2 ± V4e 2 + (Jf 2 - J*)2 ^ (Jf 2 + Jf 2 ) 2 . 

Examples of two-qubit measurements. — Using Eq. l|5}, 
we can obtain the two-qubit operations by choosing system- 
specific parameters. For example, the two-qubit operation 
with fixed e z for the Heisenberg model, XX Z model, and 
the XY model can be obtained from Eq. l|5) by setting pa- 
rameters Jf nn = J^n — Jmn> ^mn = ^mn 7^ ^mn anC ' 

Jmn = J mm J mn = 0> respectively. If all one-qubit param- 
eters are switchable, then the two-qubit operation can be ob- 
tained from Eq. (|5j by only setting e z = 0. Other effective 
spin quantum computing models presented up to now can be 
reduced by single-qubit operations to Eq. Q. For instance, 
i) the two-qubit operations of the superconducting charge 
qubit 1 12] can be reduced to Eq. (|5} with Jf 2 = Jf 2 = by a 
conjugation-by-(7r/4)(cri y + <T2y) 12511 on the Hamiltonian; ii) 
the two-qubit operation for the models in Refs. 11711 and 112 611 
can be reduced to Eq. (|5Jl with Jf 2 = Jf 2 = and model 11211 
by the conjugation-by-(7r/4)(crij / + azy) and conjugation-by- 
(7r/4)(<7i x +a2x) on the Hamiltonian of the system. Combin- 
ing the basic two-qubit operations U12 (t) with the single-qubit 
operations, we can obtain any desired two-qubit operation by 
choosing the evolution time t and the system-specific param- 
eters {e Q , Jf 2 : a — x,y, z}. 

Now, let us consider the XY exchange coupling system 
with switchable one-qubit parameters as an example to an- 
swer how to obtain the expectation values of two-qubit mea- 
surements. If we want to obtain, for example, the expectation 
value r zy of a% z <g> (72y in such system, we can first switch off 
all the one-qubit parameters £1 a , then let the two-qubit system 



evolve during the time r = far/8 Jf 2 with the evolution oper- 
ator £/i2(V), then switch off the exchange coupling and only 
make the first qubit have a ir/2 rotation Y\ = exp{— iir<Ti y /4} 
around the y axis, that is 



P 



U 12 (t) 



U 12 (t) p u\ 2 (t) ^ Y 1 U 12 (t) P u1 2 (t)Y? = p. 



(6) 

Afterwards we make the measurement (|l)(l|)i on the ro- 
tated state p obtaining the probability p = Tr[p (|1) (l|)i] = 
(\[2 + r x o + r zx )/2y2, corresponding to the equivalent two- 
qubit measurement — a\ x — <j\ z <g> a 2x . Because r X Q has been 
obtained by the equivalent single-qubit measurement <j\ x , r zx 
is completely determined by the above result. Eight other 
values of equivalent two-qubit measurements for this pair 
can also be obtained by projecting (|l)(l|)i on the measured 
state with the quantum operations summarized in Table U 
Each measured value is related to the expectation values of 
a single-qubit and a two-qubit measurements. For a two-qubit 
state in this system, the above 9 two-qubit and 6 single-qubit 
measurements are enough to obtain 15 unknown parameters 
fhh (h,h = 0,x,y,z) where h,l2 are not simultaneously 
taken as zero. The 16 matrix elements of the two-qubit state 
are obtained by combining the 15 parameters with the 
normalization condition and finally the two-qubit state can be 
completely reconstructed. 

The implementation of equivalent two-qubit measurements 
with a well-chosen two-qubit operation for a pair of cou- 
pled two-qubit system plays a significant role in the recon- 
struction of a state. For the XY and Heisenberg models 
with switchable one-qubit parameters, the equivalent mea- 
surements \f2W^ a\ z W and 2W^ai z W, to obtain the expec- 
tation values of 9 two-qubit measurements, are summarized 
in table H] where the non-local two-qubit operation operators 
Ui and U2 for the XY and Heisenberg models are chosen by 
Eq. (|5} with the system specific parameters and the evolution 
time r = hir/8 Jf 2 as 



2y/2Ui 
2V2U2 



V2 + 1 ) / + f V2 - 1 ) a\ 



iul <8> a 2 



v v x 
(2-i)I-ial 



(7) 



(8) 



The reconstructions of the qubit states in these models with 
fixed e z are the same with switchable one-qubit parameters, if 
the ratios e z j ' Jf 2 = 4m/ (2n— 1) (m, n = 1, 2, • • •) are appro- 
priately chosen and the operation time is r = (ntm) / (2e z ). 

We also find that the tomographic measurement steps for 
most systems can be reduced to the same steps needed for 
the XY model. For example, i) by choosing appropriate val- 
ues of Jf 2 , Jf 2 (^12) J\2* an< ^ £ z) and operation time t for 
the XX Z model with the switchable one-qubit parameters 
(fixed e z ) such that Jf 2 r/h = 2mr, Jf 2 r/h = (2m — 1)77/8 
(J1 2 t/H = 2mr, J? 2 t/H = (2m - 1)tt/8, £ z t/H = Itt/2) 
with l,m,n = 1, 2 ■ ■ ■ , then we can obtain the same two- 
qubit operation as for the XY model and the qubit state can 
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be reconstructed by using the same steps as the XY model, 
ii) the qubit state of the superconducting charge qubit model 
can also be reconstructed by using the same steps as the XY 
model when the parameters and evolution time are appropri- 
ately chosen 12711 , in) the qubit states in the systems modelled 
in Refs fl7il and i2al can also be reconstructed by using the 
same steps used for the XY model. 

It should be emphasized that the different qubit measure- 
ments on the quantum state with fixed quantum operations 
produce different results and the quantum operations are not 
unique to obtain each expectation value of the measurement. 
In table [I] we only discuss the procedure for the first qubit 
POM, if we can make all single-qubit measurements, then 
the operation steps to obtain some of the expectation val- 
ues of the multiple-qubit measurements can be decreased. 
For example, if we can experimentally perform the second 
qubit projection ( 1 1 ) ( 1 1 ) 2 in the XY model, the expectation 
value of <7\ y (g> CT2z can be obtained by using two steps of 
operations U\2{t) and Yjj, that is U\ 2 {t)Y2 °2z Y2U\2{t) — 
— {pix + <*\y ® o'2z)/v / 2- But four steps are needed for the 
first qubit measurement. The price paid is that noise may be 
increased because the system is in contact with more probes. 

Multi-qubit measurements. — These measurements can also 
be obtained by designing step-by-step single- and two-qubit 
operations. In principle, to obtain the expectation value of an 
/-qubit (2 < I < n) measurement, we need at least I — 1 two- 
qubit operations for different pairs of the Z-qubit. For example, 
let us obtain the expectation value r zzx corresponding to the 
three-qubit measurement <j\ z (8 <J2z <8 <J3 X for the XY inter- 
action system with switchable one-qubit parameters. We can 
replace the first qubit by the third qubit in the two-qubit oper- 
ation Ui2(t) and perform an operation U2z{t) on the second 
and third qubits, followed by another operation Uv2{t) on the 
first and second qubits, followed by a 7r/2 rotation Y\ about 
the y axis for the first qubit, followed finally by the measure- 
ment (|l)(l|)i. That is, an equivalent three-qubit measure- 
ment can be obtained as 

UL (r) U\ 2 {r)Ylu Xz Y x U 12 (r) U 23 (r) 

= - T-7=<Tlx ~ -7°lz ® °2y + ~0\z ® °2z ® &3x, (9) 

2V2 4 4 

where the assumption of exchange couplings are the same for 
all qubit pairs used. The probability of measuring (|1)(1|) 1 
on the above rotated three-qubit state is p' = (2\/2 + 2r x oo + 
V^r zy o — V2r zzx )/iy/2. Then r zzx can be determined by the 
p, r x oo, and r zy o, the latter two parameters have been obtained 
by single- and two-qubit measurements. We can also obtain 
other probabilities of the equivalent three-qubit measurements 
related to the expectation values of the three-qubit measure- 
ments by projecting (|1)(1|); on the final operated state. Fora 
three-qubit state, we can solve the equations for all probabil- 
ities of equivalent one-, two-, and three-qubit measurements 
to obtain expectation values of all measurements, finally all 
matrix elements of a three-qubit state are obtained by these 
expectation values, and the state is reconstructed. Any n- 
qubit measurement can be obtained in a similar way to the 



three-qubit measurement, then the n-qubit state can also be 
reconstructed. 

Discussion. — In summary, we have proposed a scheme 
for tomographic reconstruction of qubit states for a class of 
promising solid-state quantum computing models. We find 
that elemental logic gates, such as CNOT gate, control phase 
gate, etc. are not necessary in this process. An appropriate 
non-local two-qubit operation is enough to realize this pur- 
pose. The generalization of the above discussion to the spin 
measurement 1 16, 17] is straightforward because of the equiv- 
alence between |1)(1| = h(cro ~ °>) an d G x (Py) by a 7r/2 
(—it/2) rotation about the y (x) axes. Using present technol- 
ogy, our proposal is experimentally feasible in these solid state 
qubit systems. Ideally, the reconstructed qubit state p should 
satisfy the properties of the normalization, positivity, and Her- 
miticity. However, we always deal with a limited amount of 
experimental data, which are also affected by noise and im- 
perfect quantum measurements. To overcome problems due 
to unavoidable statistical errors, we can use the maximum- 
likelihood estimation of the density matrix 1 28] to obtain a 
more accurate reconstructed qubit. 
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